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abstract 

We study the three-dimensional (3D) compact U(l) lattice gauge theory coupled with iV-flavor 
Higgs fields by means of the Monte Carlo simulations. This model is relevant to multi-component 
superconductors, antifcrromagnetic spin systems in easy plane, inflational cosmology etc. It is known 
that there is no phase transition in the N = 1 model. For N = 2, we found that the system has a 
second-order phase transition line £1(02) in the C2(gauge coupling)— c\ (Higgs coupling) plane, which 
separates the confinement phase and the Higgs phase. Numerical results suggest that the phase 
transition belongs to the universality class of the 3D XY model as the previous works by Babaev et 
al. and Smiseth et al. suggested. For N = 3, we found that there exists a critical line similar to that 
in the N = 2 model, but the critical line is separated into two parts; one for C2 < C2tc = 2.4 ± 0.1 
with first-order transitions, and the other for C2tc < ci with second-order transitions, indicating the 
existence of a tricritical point. We verified that similar phase diagram appears for the N = 4 and 
N = 5 systems. We also studied the case of anistropic Higgs coupling in the N — 3 model and 
found that there appear two second-order phase transitions or a single second-order transition and 
a crossover depending on the values of the anisotropic Higgs couplings. This result indicates that 
an "enhancement" of phase transition occurs when multiple phase transitions coincide at a certain 
point in the parameter space. 
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1 Introduction 



There are many interesting physical systems involving multi-component (iV-componcnt) matter 
fields. Sometimes they are associated with exact or approximate symmetries like "flavor" symmetry. 
In some cases, the large- N analysis or the 1/N expansion[l] for an iV-flavor system is applicable and 
it gives us useful information that cannot be obtained by the ordinary perturbative calculations. 
But the properties of the large- N systems may differ from those at medium values of N that one 
actually wants to know. Study of the A^-dependence of various systems is certainly interesting, but 
has not been examined well. 

Among these "flavor" physics, the effect of matter fields upon gauge dynamics is of quite general 
interest in quantum chromodynamics, strongly correlated electron systems, quantum spins, etc. [2, 
3, 4, 5]. In the present paper, we shall study the three-dimensional (3D) U(l) gauge theory with 
multi-component Higgs fields <p a (x) = \<p a (x)\ exp(iip a (x)) (a = l,---,N) with fixed amplitudes 
|0 a (x)| = 1. This model is of general interest, and knowledge of its phase structure, order of its 
phase transitions, etc. may be useful to get better understanding of various physical systems. These 
systems include the following: 

N -component superconductor: Babaev[6] argued that under a high pressure and at low temper- 
atures hydrogen gas may become a liquid and exhibits a transition from a superfluid to a super- 
conductor. There are two order parameters; <j) e for electron pairs and 4> p for proton pairs. They 
may be treated as two complex Higgs fields (N = 2). In the superconducting phase, both <j) e and 
(j> p develop an off-diagonal long-range order, while in the superfluid phase, only the neutral order 
survives; lim| a .|_ >oo (<^ e (a;)^p(0)) ^ 0. 

p-wave superconductivity of cold Fermi gas: Each fermion pair in a p-wave superconductor has 
angular momentum J — I and the order parameter has three components, J z = —1,0,1. They 
are regarded as three Higgs fields (N — 3). As the strength of attractive force between fermions is 
increased, a crossover from a superconductor of the BCS type to the type of Bose-Einstein conden- 
sation is expected to take place [7]. 

Phase transition of 2D antiferromagnetic(AF) spin models: In the s = 1/2 AF spin models, 
a phase transition occurs from the Neel state to the valence-bond solid state as parameters are 
varied. Senthil et al.[8] argued that the effective theory describing this transition take a form of 
U(l) gauge theory of spinon (CP 1 ) field z a (x) (\zi\ 2 + \z2\ 2 = 1) with an additional Berry phase 
and a "deconfined phase" of spinons appears at the critical point. If one consider the easy-plane 
limit (S z — 0), \zi\ 2 — 1 1 2 = 1/2 and the CP 1 bosons are expressed by two Higgs fields as 
z a = cxp(i</? a )/v / 2 (N = 2). This system is studied by Nogueira et al. by a renormalization group 
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analysis [9] and it is clarified that the above easy-plane limit has only a first-order phase transition. 

Similar limit to the above may be taken to relate the superconductivity of ultracold fermionic 
atoms with spin J to the U(l) gauge model with N Higgs fields. J. Zhao, et al. considered the 
SU(N) Hubbard model to describe the superconductivity of fermionic atoms, which has a N = 
2 J + 1-component order parameter [10]. At large repulsion U and at the filling factor n = l/N, the 
model becomes the U(l) gauge model with CP N ~ X spins. A CP^ 1 variable is parametrized as 
z a = pa exp(i<y2 a ) with Y^l=i Pa = 1- I n the symmetric limit, which is the easy-plane limit for N = 2, 
p\ = l/N and z a becomes a Higgs field. 

Effects of doped fermionic holes (holons) to these AF spins were also studied extensively. The 
effective theory obtained by integrating out holon variables may be a U(l) gauge theory with N = 2 
Higgs fields (with nonlocal gauge interactions). Kaul et al.[ll] predicts that such a system exhibits 
a second-order transition, while numerical simulations of Kuklov et al.[12] exhibit a weak first-order 
transition. This point should be clarified in future study. 

Inflational cosmology: In the inflational cosmology [13], a set of Higgs fields is introduced to 
describe a phase transition and inflation in early universe. Plural Higgs fields are necessary in a 
realistic model [14]. 

In the rest of the present paper, we shall study the multi-Higgs lattice models by Monte Carlo(MC) 
simulations. We consider the simplest form of the model, i.e., the 3D compact lattice gauge theory 
without Berry's phase. We introduce the Higgs fields 4> xa on the site x of the cubic lattice and treat 
them in the London limit, \(j) xa \ = 1, 4>xa — exp(iip xa ). We also put the compact U(l) gauge field 
U Xfl = exp(i6 XIJ/ ) on the link (x, x + fi). p(= 1, 2, 3) is the direction index (we use them also as the 
unit vectors). The action S consists of the Higgs coupling with its coefficients c la (a = 1, . . . , N) 
and the plaquette term with its coefficient c 2 as 



A couple of models close to Eqs.(l.l) and (1.2) have been investigated. Smiseth et al.[15] studied 
the noncompact U(l) Higgs models. A duality transformation maps the charged sector into the 
inverted XY spin model. Thus they predicted that the system exhibits a single inverted XY tran- 
sition and N — 1 XY transitions. Their numerical study confirmed this prediction for N — 2. For 



1 

S = 2EE( c ^U,a^- + H - c -)+f E (ULUI + ^U W U^ + R.c). (1.1) 




The partition function Z of the model is given by 




(1.2) 
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N = 2, Kragset et al.[16] studied the effect of Berry's phase term in the compact Higgs model. They 
reported that Berry's phase term suppresses monopoles (instantons) and changes the second-order 
phase transitions to first-order ones. 

The phase structure of the present system (1.1) can be studied by the following consideration 
developed by Smiseth et al.[15, 17]. That is, among N phases ip a (x) of the Higgs fields, the sum 
(p + = J2 a ip a couples to the gauge field and describes charged excitations, whereas the remaining 
N — 1 independent linear combinations (pi{i = 1, • • • , N — 1) describe neutral excitations. The latter 
N — 1 modes may be regarded as a set of N — 1 XY spin models. As the N = 1 compact U(l) 
Higgs model stays always in the confinement phase[18], we expect N — 1 second-order transitions of 
the type of the XY model. The above discussion is useful to get an intuitive picture of the phase 
structure. In the present paper, we shall study the system (1.1) by means of MC simulation and 
verify the above conclusion. 

The rest of the present paper is organized as follows. In Sec. 2, we shall study the N = 2 multi- 
Higgs model and report the results of the MC simulations, in particular, its phase diagram in the 
Ci — c 2 plane. We measure the internal energy, specific heat and instanton density. In Sec. 3, we shall 
investigate the N = 3 cases. By varying the ratios of the three Higgs couplings, the model exhibits 
interesting phase structure. In Sec. 4, we present the result of the cases of N — 4 and N — 5. Section 
5 is devoted for conclusion. 

2 Two-flavor Higgs model (TV = 2) 

For the MC simulations, we used the standard Metropolis algorithm [19]. We consider the cubic 
lattice with the periodic boundary condition and its size L 3 up to L = 36. The typical statistics 
used was 6 • 10 5 MC steps and the averages and errors were estimated over 20 samples. 

We first study the N = 2 case with symmetric couplings Cn = C12 = c\. We measured the 
internal energy U = — (S)/L 3 and the specific heat C = ((S— (S)) 2 ) / L 3 in order to obtain the phase 
diagram and determine the order of phase transitions. 

In Fig. 1(a), we show C at c 2 = 0.4 as a function of c\ for L = 22, 28, 32. The peak of C develops 
as the system size is increased. The results indicate that a second-order phase transition occurs at 
ci ~ 0.91. In fact, we applied the finite-size-scaling (FSS) hypothesis to C of Fig. 1(a) in the form of 

C(c u L) = L^r ) (L 1 ^e), (2.1) 

where e = (ci — cio^/ci^ and Cloc is the critical coupling at L — > oo. We determined v = 
0.65 ± 0.02, a = 0.16 ± 0.02, and c loc = 0.909 ± 0.010 with the scaling function T](x) plotted in 
Fig. 1(b). This result supports the FSS (2.1). 
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Figure 1: (a) System-size dependence of specific ] 
(b) Scaling function rj(x) of Eq.(2.1) for Fig. 1(a). 
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; C of L = 22, 28, 32 for TV = 2 along c 2 = 0.4. 



The above results for N — 2 are consistent with the prediction given in the introduction. The 
sum (p x+ = (p x i + (p X 2 couples with the compact gauge field and generates no phase transition [18], 
while the difference Cp x _ = ip x \ — (p x2 behaves like the angle variable in the 3D XY model. The 3D 
XY model has a second-order phase transition with the critical exponent v = 0.666. ..[20]. Our value 
of v obtained above is consistent with this value. In fact, the same result was previously obtained 
in Ref. [16], which studied the N = 2 model on the specific line c\ — c 2 in the c 2 — c\ plane by means 
of the MC simulations of large system sizes 1 . 

It is instructive to see the behavior of the instanton density p in order to study the gauge 
dynamics at the phase transition point. We employ the definition of p in the 3D U(l) compact 
lattice gauge theory given by DeGrand and Toussaint[22]. p in Fig. 2 decreases very rapidly near the 
phase transition point at c\ ~ 0.91. This indicates that a "crossover" from dense to dilute instanton 
"phases" is accompanied with the phase transition. In other words, the observed phase transition 
can be interpreted as a confinement (small ci)-Higgs (large ci) phase transition. 

The above conclusion is supported by the following consideration. An effective gauge model 
S c g(U) is obtained by integrating out the Higgs fields <p xa in Z of Eq.(1.2), 

Z = y"[dt/]exp[5 cff (C/)], 

1 The N = 2 model with a noncompact gauge action was studied by Motrunich et al.[21]. Phase stucture of that 
mode was also clarified by the paper by Kragset et al. [16]. 
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Figure 2: Instanton density p for N — 2 at c 2 = 0.4 as a function of c\. System size L = 16. 



For small c\, the above integration over the Higgs fields can be performed by using the hopping 
expansion in powers of c\. The resultant effective action S e s(U) contains nonlocal interaction terms 
of the gauge field U xfl . Recently we studied models of U(l) gauge field in 3D, which contain nonlocal 
interactions [3]. We found that the nonlocal terms give dominant effect on the gauge dynamics and 
a confincmcnt-deconfinement phase transition takes place as their coefficients are getting large. 

Knowledge of the phase structure of the CP N ^ model is also useful to identify the phases in 
the present model. As explained in the introduction, the present model is the easy-plane limit of 
the CP 1 model. The CP^ -1 model in 3D was studied both analytically by means of the 1/N 
expansion [23] and numerically by defining the model on the lattice [24, 4] . These studies show that 
the spontaneous breaking of the internal SU(N) symmetry accompanys the phase transition to the 
Higgs phase of the gauge dynamics. Similarly in the present model, the observed second-order phase 
transition corresponds to the transition from the global U(l) symmetric phase (confinement phase) 
to the phase of the the spontaneous breaking of the global U(l) symmetry (Higgs phase). 

In Fig. 3, we present the phase diagram for N — 2 in the C2-C1 plane. There exists a second-order 
phase transition line separating the confinement and the Higgs phases. There also exists a crossover 
line similar to that in the 3D N — 1 U(l) Higgs model[18]. 
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Figure 3: Phase diagram for N = 2. There are two phases, confinement and Higgs, separated by a 
second-order phase transition line. There also exists a crossover line in the confinement phase sepa- 
rating dense and dilute instanton-density regions. Location of the phase transition line is determined 
from the data of L = 24. 

3 Three-flavor Higgs model (N = 3) 

3.1 Symmetric case c\ a = C\ 

Let us turn to the N — 3 case. Among many possibilities of three ci a 's, we first consider the 
symmetric case en = C12 = C13 = c\. One may expect that there are two (N — 1 = 2) second-order 
transitions that may coincide at a certain critical point. Studying the N = 3 case is interesting from 
a general viewpoint of the critical phenomena, i.e., whether coincidence of multiple phase transitions 
changes the order of the transition. We studied various points in the ci — c\ plane and found that 
the order of transition changes as C2 varies. 

In Fig. 4, we show U and p along C2 = 1.5 as a function of c\. Both quantities show hysteresis 
loops, which are signals of a first-order phase transition. In Fig. 5, we present C at ci = 3.0. The 
peak of C at around c\ ~ 0.48 develops as L is increased, whereas U shows no discontinuity and 
hysteresis. Therefore, we conclude that the phase transition at (c2,ci) ~ (3.0,0.48) is of second 
order. 

In order to locate the tricritical point from the first to second phase transitions, we studied 
the region ci — 2.0 ~ 2.5 in detail. In Fig. 6, we present the internal energy U and the specific 
heat C at C2 = 2.0, 2.2, 2.3, 2.4, 2.5 for the system size L = 32. U at C2 = 2.0 shows a hysteresis 
at ci ~ 0.506 ~ 0.508. As C2 increases, the hysteresis becomes milder and at C2 = 2.5 it almost 
disppears. Also, as C2 increases, the region and height of the peak of C measured from the smooth 
background become reduced. These behavior, together with the size dependence of U and C, suggest 
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Figure 4: (a) Internal energy U and (b) instanton density p for N = 3 at c% = 1.5 and L = 16. 
Both exhibit hysteresis loops in the path where c\ is first increased and then decreased by the step 
Aci = 0.0005, indicating a first-order phase transition at c\ ~ 0.551. 
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Figure 5: (a) Specific heat for N = 3 at C2 = 3.0. (b) Close-up view near the peak. The peak 
develops as L increases. 
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Figure 6: (a) Internal Energy U and (b) Specific heat C for N — 3 and L — 32 at C2 = 2.0 ~ 2.5. 
(Values of ci are indicated near each curve.) As C2 increases, the changes of U become milder and 
the region and size of the peak of C become reduced. 
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that the change in the order of phase transition from the first-order one to the second-order one 
takes place in this region of c 2 . To support this point, we also studied the distribution p(E) exp(— E) 
of the internal-energy, which is defined as 



[dU][d<j)}exp(S) 



dE J[dU][d(f>}exp{S)S{S + E) = J dE p{E) exp(-E). 



(3-1) 



If the phase transition is first order, p{E) exp(— E) exhibits a double peak structure near the critical 
point, while a second-order transition exhibits a single peak structure. In Fig. 7, we present the 
distribution of p(E) exp(-E) for c 2 = 2.0 ~ 2.5 and L = 32. 
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Figure 7: Distribution of the energy, p(E) exp(— E) of Eq.(3.1), for N = 3 and L = 32. (a) 
(c 2 , Cl ) = (2.0,0.507); (b) (c 2 , Cl ) = (2.3,0.495); (c) (c 2 , Cl ) = (2.4,0.493); (d) (c 2 , Cl ) = (2.5,0.492). 
The signal of double-peak structure, which is shown clearly in Fig. (a), becomes weaker as c 2 increases, 
and disappears in Fig.(d). 
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The value of c\ is chosen near the peak location of C given in Fig. 6(b). Apparently, Fig. 7(a), the 
distribution at C2 = 2.0, fits better by a double-peak (e.g., double Gaussian) distribution rather than 
by a single-peak one. The fact that two peaks here have different weights mainly reflects that it is 
slightly away from the critical point of c\. Fig. 7(c) for C2 = 2.4 has a single peak at E ~ —7.73, 
but not symmetric around this peak, which shows a remnant of the second peak at lower U region. 
Fig. 7(d) for C2 = 2.5 shows a clean single-peak distribution. From these observations we determine 
that the tricritical point is located at C2 = C2t c — 2.4 ± 0.1. 

In Fig. 8 (a), we present the phase diagram of the symmetric case for N — 3, where the order of 
transition between the confinement and Higgs phases changes from first (c 2 < Cat c ) to second order 
(c2tc < C2). In Fig. 8(b) we present C along c\ = 0.2, which shows a smooth nondcvcloping peak. As 
shown in Fig. 9, the instanton density p decreases smoothly around this peak. These results indicate 
a crossover at ci ~ 1.5. 




Figure 8: (a) Phase diagram for the N = 3 symmetric case. The phase transitions are first order 
in the region C2 < C2tc — 2.4, whereas they are second order in the region C2 > C2tc- There exists 
a tricritical point at around (c2,ci) ~ (2.4,0.49). Crosses near C2 = 1.5 line show crossovers, (b) 
Specific heat for N = 3 at ci = 0.2. It has a system-size independent smooth peak at which a 
crossover takes place. 



3.2 Asymmetric case 

Then it becomes interesting to consider asymmetric cases, e.g., cu 7^ C12 = C13. This case is closely 
related to a doped AF magnet. 4>2 and ^3 correspond there to the CP 1 spinon field in the deep easy- 
plane limit, whereas <pi corresponds to doped holes (although they are fermionic). This case is also 
relevant to cosmology because the order of Higgs phase transition in the early universe is important 
in the inflational cosmology. Furthermore, one may naively expect that once a phase transition to 
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Figure 9: Instanton density p in the N — 3 symmetric case for c\ = 0.2 as a function of C2- System 
size L = 16. 



the Higgs phase occurs at certain temperature T, no further phase transitions take place at lower 
T's even if the gauge field couples with other Higgs bosons. However, our investigation below will 
show that this is not the case. 

Below we shall consider the two cases, (i) Cxi — C13 > en and (ii) C12 = C13 < c\\. Let us first 
consider the case (i) cyi = C13 = 2cn, which we call the c\ = (1,2,2) model, and focus on the 
case C2 = 1.0. As shown in Fig. 10(a), C exhibits two peaks at en ~ 0.35 and 0.52. Figs. 10(b), (c) 
present the detailed behavior of C near these peaks, which show that the both peaks develop as L 
is increased. We conclude that both of these peaks show second-order transitions. This result is 
interpreted as the first-order phase transition in the symmetric N = 3 model is decomposed into 
two second-order transitions in the c\ = (1, 2, 2) model. 

Let us turn to the opposite case (ii), C12 = C13 = 0.5cn, i.e., the c\ — (2, 1, 1) model at C2 = 1.0. 
One may expect that two second-order phase transitions appear as in the previous c\ = (1,2,2) 
model. However, the result shown in Fig. 11 indicates that there exists only one second-order phase 
transition near en ~ 1.08. The broad and smooth peak near c\\ ~ 0.85 shows no L dependence 
and we conclude that it is a crossover. This crossover is similar to that in the ordinary N = 1 
gauge-Higgs system as we shall see by the measurement of p below. 

The orders of these transitions are understood as follows: In the c\ = (1,2,2) model, as we 
increase en, the two modes 4> xa (a — 2,3) with larger c\ a firstly become relevant and the model 
is effectively the symmetric N = 2 model. The peak in Fig. 10(b) is interpreted as that of the 
second-order phase transition in this model. As the Higgs couplings en's are increased further, the 
fluctuations of the gauge field is negligibly small, and the effective model is the N = 1 XY model 
of 4> x i. It gives the second-order peak in Fig. 10(c). Similarly, in the c\ — (2, 1, 1) model, <fi x i firstly 
becomes relevant. The effective model is the N = 1 model, which gives the broad peak in the specific 
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Figure 10: (a) Specific heat of the c\ — (1, 2, 2) model (iV=3) at C2 = 1.0. (b,c) Close-up views of C 
near (b) cn ~ 0.35 and (c) cn ~ 0.52. 

heat C in Fig. 11 as the crossover takes place there[18]. For larger values of en's, the effective model 
is the N = 2 symmetric model of 4> X 2j4'x3 and U xfl , giving the sharp second-order peak in C in 
Fig.ll. 

In Fig. 12, we present p of the c\ = (1,2,2) and (2,1,1) models as a function of cn. p of the 
ci = (1,2,2) model decreases very rapidly at around cn ~ 0.35, which is the phase transition 
point in lower cn region. On the other hand, at the higher phase transition point, cn ~ 0.52, p 
shows no significant changes. This observation indicates that the lower-cn phase transition is the 
confinement-Higgs transition, whereas the higher-cn transition is a charge-neutral XY-type phase 
transition. 

On the other hand, p of the C\ — (2, 1, 1) model decreases rapidly at around cn ~ 0.85, where C 
exhibits a broad peak. This indicates that the crossover from the dense to dilute-instanton regions 
occurs there just like in the N = 1 case [18]. No "anomalous" behavior of p is observed at the critical 
point cn ~ 1.1, and therefore the phase transition is that of the neutral mode. 
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Figure 11: Specific heat of the c\ — (2, 1, 1) model at C2 = 1.0. 




Figure 12: Instanton density p at C2 = 1.0 in the (a) c\ = (1, 2, 2) model and (b) c\ = (2, 1, 1) model. 
System size L = 16. 

4 Symmetric model with TV = 4 and 5 

We have also studied the symmetric case for N — 4 and 5 multi-Higgs models at C2 = 0. Both 
cases show clear signals of first-order transitions at c\ ~ 0.89(iV = 4),0.84(7V = 5) as shown in 
Fig. 13. On the other hand, at ci = oo, the gauge dynamics is "frozen" to U xfl = 1 up to gauge 
transformations, so there remain iV-fold independent XY spin models, each of which exhibits a 
second-order transition at c± ~ 0.46. Thus we expect a tricritical point for general N > 2 at some 
finite ci separating first-order and second-order transitions. 

5 Conclusion 

In the present paper, we studied the U(l) multi-flavor Higgs model in 3D, which is closely related 
to various interesting physical systems. By means of the MC simulations, we clarified its phase 
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Figure 13: Internal energy U in (a) the N = 4 model and (b) N = 5 model for C2 = with the system 
size L = 24. Both exhibit hysteresis loops in the path where ci is first increased and then decreased 
by the step Aci = 0.005. This indicates a first-order phase transition for C2 = at c\ ~ 0.89 for 
N = 4 and at ci ~ 0.84 for TV = 5. 

structure and critical behavior. Let us summarize the results. For N = 2 there is a critical line 
Cx{c%) of second-order transitions in the ci — c\ plane, which distinguishes the Higgs phase (c\ > c\) 
and the confinement phase (ci < c\). This result is consistent with Kragset et al.[16]. We obtained 
the crtitical exponent of the phase transition by means of the FSS and found that the result is very 
close to that of the 3D XY model. 

For N = 3 there is a similar transition line, but the region < ci < cite — 2.4 is of first-order 
transitions while the region C2t c < C2 is of second-order transitions. We concluded that there exists 
a tricritical point. It is very interesting and also important to clarify the nature of the tricritical 
point, especially its critical exponent. This problem is under study and the result will be published 
in future . 

To study the mechanism of generation of these first-order transitions, we studied the asymmetric 
cases and found two second-order transitions [in the c\ — (1,2,2) model] or one crossover and one 
second-order phase transition [in the c\ — (2, 1, 1) model]. The former case implies that two simulta- 
neous second-order transitions strengthen the order to generate a first-order transition. Chernodub 
et al.[25] reported a similar generation of an enhanced first-order transition in a related 3D Higgs 
model with singly and doubly charged scalar fields. We stress that the above change of the order 
is dynamical because (1) It depends on the value of C2, (2) Related 3D models, the CP N ~ 1 and 
TV-flavor CP 1 gauge models, exhibit always second-order transitions (See Ref. [4]). 
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